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Abstract 

In this paper, we propose compactly supported radial basis functions for solving some well- 
known classes of astrophysics problems categorized as non-linear singular initial ordinary dif¬ 
ferential equations on a semi-infinite domain. To increase the convergence rate and to decrease 
the collocation points, we use the compactly supported radial basis function through the integral 
operations. Afterwards, some special cases of the equation are presented as test examples to 
show the reliability of the method. Then we compare the results of this work with some results 
and show that the new method is efficient and applicable. 

Keywords: Lane-Emden type equation. Compact support radial basis functions. Isothermal Gas 
sphere. White Dwarf equation. Non-linear ordinary differential equation 


1. Introduction 

In recent decades, the so-called meshless methods have been extensively used to find ap¬ 
proximate solutions of various types of linear and non-linear equations (Fasshauer, 2007) such 
as differential equations (DEs) and integral equations (lEs). Unlike the other methods which were 
used to mesh the domain of the problem, meshless methods don’t require a structured grid and 
only make use of a scattered set of collocation points regardless of the connectivity information 
between the collocation points. 

For the last years, the radial basis functions (RBFs) method was known as a powerful tool 
for the scattered data interpolation problem. One of the domain-type meshless methods is given 
in (Kansa, 1990a) in 1990, which directly collocates RBF, particularly the multiquadric (MQ) 
to find an approximate solution of linear and non-linear DEs. The RBFs can be compactly and 
globally supported. The global RBF’s are infinitely differentiable and contain a free parameter c, 
called the shape parameter (Buhman, 2000; Islam and S. Haqb, 2009; Sara, 2005). The interested 
reader is referred to the recent books and paper by Buhmann (Buhman, 2000, 2004) and Wend- 
land (Wendland, 2005) for more basic details about RBFs, compactly and globally supported and 
the convergence rate of the RBFs. In RBF method using globally RBFs, if c increases, the system 
of equations to be solved becomes ill-conditioned. Cheng et al. (Cheng et al., 2003) showed that 


* Corresponding author. Member of research group of Scientific Computing. Fax: +98 2122431653. 

Email addresses: k_parand@sbu. ac. ir (K. Parand), mohammadheinainiOyahoo. com (M. Hemami) 

Preprint submitted to Int. J. Appl. Comput. Math May 31, 2016 





when c is very large then the RBFs system error is of exponential convergence. To overcome 
the problems of globally RBF, alternative ways suggested such as domain decomposition(Lee 
and Hon, 2004), LU decomposition(Kansa and Hon, 2000), good matrix pre-conditioners(Kansa 
and Hon, 2000) change the global support RBFs with compactly support RBFs (CSRBFs)(Shen, 
2012) which are local support and have not any parameter. 

There are two basic approaches for obtaining basis functions from RBFs, namely direct ap¬ 
proach (DRBF) based on differential process (Kansa, 1990b) and indirect approach (IRBF) based 
on an integration process (Mai-Duy, 2005; Mai-Duy and Tran-Cong, 2001a,b). Both approaches 
were tested on the solution of second order DEs and the indirect approach was found to be supe¬ 
rior to the direct approach (Mai-Duy and Tran-Cong, 2001a). 

In this paper, we use the indirect CSRBF (ICSRBF) for finding the solution of Lane-Emden 
type equations and also Isothermal gas sphere. White Dwarf equation. 

This paper is arranged as following: 

In Section 2 we describe the Lane-Emden equations. In Section 3 we survey several methods 
that have been used to solve Lane-Emden type equations. In Section 4, the properties of CSRBF 
and the way to construct the ICSRBF method for this type of equations are described. In Section 
5 the proposed method is applied to some types of Lane-Emden equations, and a comparison is 
made with the existing solutions that were reported in other published works. Einally we give a 
brief conclusion in the last section. 


2. Lane-Emden type equations 

The Lane-Emden equation describes a variety of phenomena in theoretical physics and as¬ 
trophysics, including aspects of stellar structure, the thermal history of a spherical cloud of gas, 
isothermal gas spheres, and thermionic currents. 


Let P{r) denote the total pressure at a distance r from the center of spherical gas cloud. The 
total pressure is due to the usual gas pressure and a contribution from radiation: 


P 


+ 

3 


RT 

5 

V 


where g , T, R and v are respectively the radiation constant, the absolute temperature, the gas 
constant, and the specific volume(Agrawala and O’Regnab, 2007). Let M{r) be the mass within 
a sphere of radius r and G the constant of gravitation. The equilibrium equations for the config¬ 
uration are 

dP GM{r) 

dr r^ 


dM(r) 

dr 


— Ajipr^, 


where p is the density, at a distance r from the center of a spherical star. 
Eliminating M yields: 

1 d(r^dP\ 


Pressure and density p = v ^ vary with r and P = Kp^^'” where K and m are constants. 

We can insert this relation into Eq. (1) for the hydrostatic equilibrium condition and from this 
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rewrite equation to: 


K{m + 1) i_il 1 d I 2 dy\ _ ^ 

4-7tG ” \r^dr\ dr] ~ 

where A representing the central density of the star and y is dimensionless quantity that are both 
related to p through the following relation: 
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Inserting these relations into our previous relations we obtain the Lane-Emden equation (De- 
hghan and Saadatmandi, 2007): 


lH2dy\ 
dx\ dxj 


and with simplifying previous equation we have: 


y''(x)+ -y'(x)+y"’(x) = Q, 

X 


X >0, 


with the boundary conditions: 

y(0) = l, y(0) = 0. 

It has been claimed in the literature that only for m = 0, 1 and 5 the solutions of the Lane-Emden 
equation could be exact. Eor the other values of m, the Lane-Emden equation is to be integrated 
numerically. In this paper, we solve it for m = 1,1.5,2,2.5,3,4 and 5. 


3. Methods have been proposed to solve Lane-Emden type equation 

in years, many analytical and numerical methods have been used to solve Lane-Emden equa¬ 
tions, the main difficulty arises in the singularity of the equations at x = 0. Currently, most 
techniques which were used in handling the Lane-Emden type problems are based on either se¬ 
ries solutions or perturbation techniques. 

Bender et al. (Bender et al., 1989), proposed a new perturbation technique based on an artificial 
parameter 6, the method is often called d-method. 

Mendelzweig and Tabakin (Mendelzweig and Tabakin, 2001) used quasi-linearization ap¬ 
proach to solve the standard Lane-Emden equation. This method approximates the solution of a 
non-linear differential equation by treating the non-linear terms as a perturbation about the linear 
ones, and unlike perturbation theories is not based on the existence of some small parameter. 
He showed that the quasi-linearization method gives excellent results when applied to different 
no-nlinear ordinary differential equations in physics, such as the Blasius, Duffing, Lane-Emden 
and Thomas-Eermi equations. 

Shawagfeh (Shawagfeh, 1993) applied a non-perturbative approximate analytical solution for 
the Lane-Emden type equation using the Adomian decomposition method. His solution was in 
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the form of a power series. He used Fade approximants method to accelerate the convergence of 
the power series. 

In (Wazwaz, 2001), Wazwaz employed the Adomian decomposition method with an alter¬ 
nate framework designed to overcome the difficulty of the singular point. It was applied to the 
differential equations of Lane-Emden type. In future (Wazwaz, 2006) he used the modified de¬ 
composition method for solving analytical treatment of non-linear differential equations such as 
Lane-Emden equations. The modified method accelerates the rapid convergence of the series so¬ 
lutions, dramatically reduces the size of work, and provides the solution by using few iterations 
only without any need to the so-called Adomian polynomials. 

Liao (Liao, 2003) provided a reliable, easy-to-use analytical algorithm for Lane-Emden type 
equations. This algorithm logically contained the well-known Adomian decompositions method. 
Different from all other analytical techniques, this algorithm itself provides us with a convenient 
way to adjust convergence regions even without Fade technique. 

He (He, 2003) employed Ritz’s method to obtain an analytical solution of the problem. By 
the semi-inverse method, a variational principle is obtained for the Lane-Emden equation, which 
he gave much numerical convenience when applied to finite element methods or the Ritz method. 

Farand et al. (Farand and Khaleqi, 2016; Farand and Razzaghi, 2004a,b) presented some 
numerical techniques to solve higher ordinary differential equations such as Lane-Emden. Their 
approach was based on a rational Chebyshev and rational Legendre tau method. They presented 
the derivative and product operational matrices of rational Chebyshev and rational Legendre 
functions. 

These matrices together with the tau method were utilized to reduce the solution of these 
physical problems to the solutions of systems of algebraic equations. Also, in some paper’s 
(Farand et al., 201 la,b, 2010a, 2013b), Farand et al. applied Hermite function collocation method 
(HLC) , Bessel function collocation method and meshless collocation method based on Radial 
basis function (RBLs) to solving the Lane-Emden type equations. 

Ramos (Ramos, 2003, 2005, 2007, 2008) solved Lane-Emden equation through different 
methods. In (Ramos, 2005) he presented linearzation methods for singular initial value problems 
in second order ordinary differential equation such as Lane-Emden. These methods result in 
linear constant-coefficient ordinary differential equations which can be integrated analytically, 
thus they yield piecewise analytical solutions and globally smooth solutions (Ramos, 2007). 
Later, he developed piecewise-adaptive decomposition methods for the solution of non-linear 
ordinary differential equations. Fiecewise-decomposition methods provide series solutions in 
intervals which are subject to continuity conditions at the end points of each interval and their 
adoption is based on the use of either a fixed number of approximants and a variable step size, 
a variable number of approximants and a fixed step size or a variable number of approximants 
and a variable step size. In (Ramos, 2008), series solutions of the Lane-Emden equation, based 
on either a volterra integral equation formulation or the expansion of the dependent variable 
in the original ordinary differential equation are presented and compared with series solutions 
obtained by means of integral or differential equations based on a transformation of the dependent 
variables. 

Yousefi (Yousefi, 2006) used integral operator and converted Lane-Emden equations to in¬ 
tegral equations and then applied Legendre Wavelet approximations. In this work properties of 
Legendre wavelet together with the Gaussian integration method were utilized to reduce the inte¬ 
gral equations to the solution of algebraic equations. By his method, the equation was formulated 
on [0, 1]. 

Chowdhury and Hashim (Chowdhury and Hashim, 2009) obtained analytical solutions of the 
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generalized Emden-Fowler type equations in the second order ordinary differential equations by 
homotopy-perturbation method (HPM). This method is a coupling of the perturbation method 
and the homotopy method. The main feature of the HPM (Dehghan and Shaken, 2008b) is that it 
deforms a difficult problem into a set of problems which are easier to solve. HPM yields solution 
in convergent series forms with easily computable terms. 

Aslanov (Aslanov, 2008) constructed a recurrence relation for the components of the ap¬ 
proximate solution and investigated the convergence conditions for the Emden-Fowler type of 
equations. He improved the previous results on the convergence radius of the series solution. 

Dehghan and Shakeri (Dehghan and Shaketi, 2008a) investigated Lane-Emden equation us¬ 
ing the variational iteration method and showed the efficiency and applicability of their proce¬ 
dure for solving this equation. Their technique does not require any discretization, linearization 
or small perturbations and therefore reduces the volume of computations. 

Bataineh et al. (Bataineh et al., 2009) obtained analytical solutions of singular initial value 
problems (IVPs) of the Emden-Fowler type by the homotopy analysis method (HAM). Their 
solutions contained an auxiliary parameter which provided a convenient way of controlling the 
convergence region of the series solutions. It was shown that the solutions obtained by the 
Adomian decomposition method (ADM) and the homotopy-perturbation method (HPM) are only 
special cases of the HAM solutions. 

singh et al. (Singh et al., 2009) used the modified Homotopy analysis method for solving the 
Lane-Emden-equation and White Dwarf equation. 

Adibi and Rismani in (Adibi and Rismani, 2010) proposed the approximate solutions of sin¬ 
gular IVPs of the Lane-Emden type in second-order ordinary differential equations by improved 
Legendre-spectral method. The Legendre-Gauss point used as collocation nodes and Lagrange 
interpolation is employed in the Volterra term. 

Karimi vanani and Aminataei (Vanani and Aminataei, 2010) provide a numerical method 
which produces an approximate polynomial solution for solving Lane-Emden equations as sin¬ 
gular initial value problem. They are first, used as an integral operator and convert Lane-Emden 
equations into integral equations, then convert the acquired integral equation into a power series 
and finally, they transforming the power series into pade series form. 

Kaur et al. (Kaur et al., 2013), obtained the Haar wavelet approximate solutions for the gener¬ 
alized Lane-Emden equations. This method was based on the quasi-linearization approximation 
and replacement of an unknown function through a truncated series of Haar wavelet series of the 
function. 

Other researchers try to solve the Lane-Emden type equations with several methods. For ex¬ 
ample, Yildirim and Ozi§ (Yildinm and Ozi§, 2007, 2009) by using HPM and VIM methods, 
Benko et al. (Benko et al., 2008) by using Nystrom method, Iqbal and javad (Iqbal and Javad, 
2011) by using Optimal EL4M, Boubaker and Van Gorder (Boubaker and Van-Gorder, 2012) 
by using Boubaker polynomials expansion scheme, Da§cioglu and Yaslan (Akyiiz-Dagcioglu 
and (^erdik Yaslan, 2011) by using Chebyshev collocation method, Yiizbagi (Yiizbagi, 2011; 
Yiizbagi and Sezer, 2013) by using Bessel matrix and improved Bessel collocation method, Boyd 
(Boyd, 201 1) by using Chebyshev spectral method, Bharwy and Alofi (Bharwy and Alofi, 2012) 
by using Jacobi-Gauss collocation method, Pandey et al. (Pandey and Kumar, 2012; Pandey 
et al., 2012) by using Legendre and Brenstein operation matrix, Rismani and monfared (Ris¬ 
mani and Monfared, 2012) by using Modified Legendre spectral method, Nazari-Golshan et al. 
(Nazari-Golshan et al., 2013) by using Homotopy perturbation with Fourier transform, Doha et 
al. (Doha et al., 2013) by using second kind Chebyshev operation matrix algorithm, Carunto 
and bota (Caruntu and Bota, 2013) by using Squared reminder minimization method. Mall and 
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Chakaraverty (Mall and Chakraverty, 2014) by using Chebyshev Neural Network based model, 
Giirbiiz and sezer (Giirbiiz and Sezer, 2014) by using Laguerre polynomial, Kazemi-Nasab et 
al. (Kazemi-Nasab et al., 2015) by using Chebyshev wavelet finite difference method, Hosseini 
and Abbasbandy (Hosseini and Abbasbandy, 2015) by using combination of spectral method and 
ADM method and Azarnavid et al. (Azarnavid et al., 2015) by using Picard-Reproducing Kernel 
Hilbert Space Method 

4. ICSRBF method 

4.1. CSRBF 

Many problems in science and engineering arise in infinite and semi-infinite domains. Dif¬ 
ferent numerical methods have been proposed for solving problems on various domains such as 
FEM(Bu et al., 2015; Choi and Kweon, 2016), FDM(Bu et al., 2015; Noye and Dehghan, 1999), 
Spectral(Parand et al., 2013a, 2010b, 2016, 2013b; Rad et al., 2014) methods and meshfree meth- 
ods(Abbasbandy et al., 2014; Dehghan and Shokri, 2009a; Rad et al., 2012; Shokri and Dehghan, 
2010). The use of the RBF is one of the popular meshfree method for solving the differential 
equations (Dehghan and Shokri, 2008, 2009c; Rad et al., 2015a,b). For many years the global 
radial basis functions such as Gaussian, Multi quadric. Thin plate spline. Inverse multiqudric 
and others were used (Dehghan and Shokri, 2009b; Parand and Rad, 2013; Rashidi et al., 2014) 
to solve different DFs and interpolation. These functions are globally supported and generate a 
system of equations with ill-condition full matrix. To convert the ill-condition matrix to a well- 
condition matrix, CSRBFs can be used instead of global RBFs. CSRBFs can convert the global 
scheme into a local one with banded matrices. Which makes the RBF method more feasible for 
solving large-scale problem (Wong et al., 2002). 

Wendland’s functions 

The most popular family of CSRBF are Wendland functions. These functions were intro¬ 
duced by Holger Wendland in 1995 (Wendland, 1995). he started with the truncated power 
function (pi{r) = (1 - which is strictly positive definite and radial on for 1 > [f J + 1 , and 
then he walks through dimensions by repeatedly applying the monte operator (I). 

Definition 1 (Fasshauer, 2007) with ffr) = (1 - he defines 

‘Ps.k - l’^‘Pi^\+k+\, (2) 

it turns out that the functions 0s ^ are all supported on [0,1]. 

Theorem 1 (Fasshauer, 2007) The functions fs.k are strictly positive definite (SPD) and radial 
on K* and are of the form 

(ps,k{r) re [0,1], 

10 r > 1, 

with a univariate polynomial ps,/t of degree L|J + 3k + 1. Moreover, fs.k £ are unique 

up to a constant factor, and the polynomial degree is minimal for given space dimension s and 
smoothness 2k (Fasshauer, 2007). Wendland gave recursive formulas for the functions 0s,/t for 
all s, k. 
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Theorem 2 (Fasshauer, 2007) The functions k - 0, 1, 2, 3, have form 
4>s,o ^ {I-r)[, 

cps,! = (1 - + 4/ + 3)r^ + (31 + 6)r + 3], 

<^.,3 = (1 - ryf[(P + 9fi + 231 + 15)r^ + (6f + 361 + 45)r^ 
+(15Z + 45)r+ 15], 


where I = and the symbol = denotes equality up to a multiplicative positive constant. 

The case k - 0 directly follows from the definition 1. application of the definition 1 for the case 
k - 1 yields 


4>s,i(r) 



(m(r) 


f(l - ty^dt — 


Jr 

r 


f(l - t^dt 


1 

(/ + 1 )(/ + 2 ) 


(i_^y+i[(/+iy+i]. 


where the compact support of (pi reduces the improper integral to a definite integral which can 
be evaluated using integration by parts. The other two cases are obtained similarly by repeated 
application of /.(Fasshauer, 2007) We showed most of the wendland functions in table 1 . 


Table 1: Wendland’s compactly supported radial function for various choices of k and s=3. 


^s,k 

smoothness 

SPD 

<PxQ(r) = (1 - r )2 

C“ 

m3 

(pi,\(r) = (1 -r)\(4r+ 1 ) 

C 2 

m3 

03,2(r) = (1 - r)6(35r2 + 18r + 3) 

& 

m3 

(Pxz(r) = (1 - r)l(32p + 25r^ + 8 r + 1 ) 

c® 

m3 

03 _ 4 (r) = (1 - r)'f(429r^ + 450r3 + 210r2 + 50r + 5) 


m3 

03 5 (r) = (1 - r)f (2048r5 + 2697/ + 1644/ + 566/ + 108r + 9) 

CW 

m3 


4.2. Interpolation by CSRBFs 

To interpolate or approximate one dimentional function we can represent it by a CSRBF 
as 


N 

y(-^) ~ yM = ^ ^i(Pi(x) = <i>^(3c)S, 

i=l 


(pi(x) = (p( 


\\x-Xi\\ 

ra> 


), 


(i>^(x) = [01 (x), 02 (x), ■ ■ ■ , (Pn(x)], 


2 , 
1 


where 









( 3 ) 


X is the input, r^, is the local support domain and ^,s are the set of coefficients to be determined. 
By using the local support domain, we mapped the domain of problem to CSRBF local domain. 
By choosing N interpolate points (xj, j - 1,2, ■ ■ ■ , A^) in domain: 

N 

yj = = 1,2, • • • ,A^). 

i=l 

To summarize the discussion on the coefficients matrix, we define 

AS = Y, (4) 


where : 


Y 

= [yi,y2, 

■■■,yN]\ 


A 

= [(l)^(^l),0^(V2),-- 



'4>i(xi) 

(I>2{X\) ■ ■ 

■ MxiV 


4>i(X2) 

(l>2(X2) ■ ■ 

■ 4>n(X2) 


-.<Pi{xn) 

4>2{Xn) ■ • 

<f>N(XN)j 


Note that ^i(xj) 
achieved. 


4 >( 



), by solving 


the system (4), the unknown coefficients will be 


4.3. ICSRBF method 

In this paper, we construct the 0,(v) in Eq. (5) by using the Wendland function with parame¬ 
ter k = 3 and I - 5. 

At first, we approximate the highest order derivative in the problem by expansion of CSRBFs: 

j 2 N 

-^yix) - y®(x) = ^ (5) 

1=1 

where (pi(x) is the CSRBF and is the local support domain and coefficients must be deter¬ 
mined. Now we define the lower order of derivatives and unknown function y(x) using Gauss- 
Fegendre integration as follows: 

^ yf’Wt - ^3'('^)U=o = ^ ^ 2 (6) 

y(x) - f - y(0) = (y„(x)), 

Jo 

where w[j] are weighted coefficients of Gauss-Legendre integration and define as follow; 


8 


(7) 






Pq{x) is Legendre polynomial of order q and q[i] is y-th root of 
To satisfy the initial conditions of problem, we have: 


f y?{t)du 

Jo 

(8) 

(yn(x))^ f (yi^\t))dt+l, 

Jo 

(9) 

then, we produce the residual function as follows: 


Res(x) = xyj^\x) + a(yi^\x)) + x(y„(x))™. 

(10) 

now, obtain A-nodes as follows: 


Xj^L(-y, 7 = 1,2,---,A, 

(11) 


where L is the last collocate node in domain and y is arbitrary parameter. For Lane-Emden type 
equations, we selected y between 1.5 and 1. 8 , because with these values we found the best ||/?es ||2 
for solution equations. 

By using Xj,j - 1,2, ...,A^ and residual function, we obtain N equation, so by solving these 
equations we obtain the unknown coefficients 

The result of this section can be summarized in the following algorithm for the IVP: 

F{x,y{x),y\x),y'\x))^Q, y(0) ^ a, y'(0)^b. (12) 

Algorithm The algorithm works in the following manner: 

(1) Choose N center points from domain [0, oo). 

(2) Approximate y"(x) as the from u'^(x) - 

(3) Obtain y'{x) by using defined integral operation I^{h{x)) = h{t)dt in the form u'j^{x) - 

fo'‘ + b. 

(4) Obtain y(x) by using defined integral operation I^(h(x)) - h{t)dt in the form m^vCx) = 
u'f^{t)dt + a. 

(4) Substitiute un(x), m^(x) and u'i^(x) into the main problem and create residual function Res(x). 

(5) Substitiute collocation points {^ 7)^1 into the Res(x). 

( 6 ) Solve the N equations with N unknown coefficients and find the numerical solution. 

Res(xj) = 0, 7 = 1,2, • • • , A. 


5. Application 

In this section, we apply ICSRBF method to solve the Lane-Emden type equations. In general 
the Lane-Emden type equation are formulated as follows: 

y"(x) +—y'(x) + p(x)q(y(x)) - h(x), ax > 0, (13) 


with initial conditions : 


y(0)=A, /(0) = B, 
9 


( 14 ) 



where a, A and B are real constants and p(x), q(y(x)) and h(x) are some given function. 
To apply the collocation method, we construct the residual function by substituting y„(x) in the 
Lane-Emden type Eq. (13): 

Res{x) = yf\x) + + p{x)q(yn{x)) - h{x). 

The equation for obtaining the coefficient arise from equalizing Res(x) to zero at N point 

( 11 ): 

Res(xj)^0, j ^0,1,2,...,N. (15) 

By solving this set of equations we obtain a approximate function y„(x). Note that these N 
equations generate a set of N nonlinear equations which can be solved by a well-known method 
such as Newton method for unknown coefficient 

5.7. Example 1 (The standard Lane-Emden equation) 


Table 2: Compaiison of the first zeros of standard Lane Emden equations between the present method and exact numerical 
value given by Horedt(Horedt, 2004). 


m 

N 

ra, 

L 

The present method 

Horedt 

Error 

0 

40 

6.5 

10 

2.44948974 

2.44948974 

O.OOe-00 

1 

40 

6.5 

10 

3.14159265 

3.14159265 

O.OOe-00 

1.5 

30 

2 

4 

3.65375388 

3.65375374 

1.40e-07 

2 

20 

4 

6 

4.35287462 

4.35287460 

2.00e-08 

2.5 

20 

4 

6 

5.35527531 

5.35527546 

1.50e-07 

3 

20 

6.5 

8 

6.89684855 

6.89684862 

7.00e-08 

4 

20 

14 

16 

14.9715808 

14.9715463 

3.45e-05 


Table 3: Comparison of y(x') of standard Lane-Emden equation between present method and exact values given by 
Horedt(Horedt, 2004), for m = 2. 


X 

ICSRBE 

Horedt 

Error 

0.00 

1.0000000 

1.0000000 

O.OOe-HOO 

0.1 

0.9983350 

0.9983349 

O.OOe-HOO 

0.5 

0.9593527 

0.9593527 

O.OOe-HOO 

1.0 

0.8486541 

0.8486541 

O.OOe-HOO 

3.0 

0.2418240 

0.2418241 

l.OOe-07 

4.0 

0.0488398 

0.0488401 

3.00e-07 

4.3 

0.0068107 

0.0068109 

2.00e-07 

4.35 

0.0003660 

0.0003660 

O.OOe-HOO 


Eor p(x) - 1, q{y{x)) - y^{x), a = 2, A = 1, h(x) - 0, and B - 0, Eq. (13) is the 
standard Lane-Emden eqution that was used to model the thermal behaviour of a spherical cloud 


10 










Table 4: Comparison of of standard Lane-Emden equation between present method and exact values given by 
Horedt(Horedt, 2004), for m = 3. 


X 

ICSRBF 

Horedt 

Error 

0.00 

1.0000000 

1.0000000 

O.OOe+00 

0.1 

0.9983358 

0.9983358 

O.OOe+00 

0.5 

0.9598391 

0.9598391 

O.OOe+00 

1.0 

0.8550575 

0.8550576 

l.OOe-07 

5.0 

0.1108198 

0.1108198 

O.OOe+00 

6.0 

0.0437379 

0.0437380 

l.OOe-07 

6.8 

0.0041677 

0.0041678 

l.OOe-07 


Table 5: Comparison of y(jr) of standard Lane-Emden equation between present method and exact values given by 
Horedt(Horedt, 2004), for m = 4. 


X 

ICSRBF 

Horedt 

Error 

0.00 

1.0000000 

1.0000000 

O.OOe+00 

0.1 

0.9983401 

0.9983367 

3.40e-06 

0.2 

0.9933921 

0.9933862 

5.90e-06 

0.5 

0.9603107 

0.9603109 

2.00e-07 

1.0 

0.8608195 

0.8608138 

5.70e-06 

5.0 

0.2359598 

0.2359227 

3.71e-06 

10 

0.05965343 

0.0596727 

1.92e-05 

14 

0.0083590 

0.0083305 

2.85e-05 

14.9 

0.00058404 

0.0005764 

7.64e-06 


Table 6: of the standard Lane-Emden equations for m = 1.5,2,2.5,3,4 respectively. 


n 

m=1.5 

m=2 

m=2.5 

m=3 

m-4 

5 

5.516e-03 

1.434e-02 

1.338e-03 

9.996e-02 

1.869e-00 

10 

7.250e-05 

1.646e-04 

7.148e-06 

1.851e-04 

6.436e-03 

15 

1.050e-06 

1.379e-05 

5.919e-09 

1.699e-08 

5.362e-03 

20 

1.731e-08 

1.869e-06 

1.268e-12 

9.786e-10 

3.049e-05 


Table 7: maximum error of the standard Lane-Emden equations for m = 0,1,5 respectively. 


n 

m=0 

m=l 

m=5 

5 

2.70199e-01 

7.87765e-02 

3.54822e-01 

10 

6.9081 le-03 

2.30291e-03 

6.61006e-03 

15 

1.15061e-03 

2.87683e-04 

9.08909e-04 

20 

2.41206e-04 

4.78916e-05 

1.71788e-05 
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of gas acting under the mutual attracting of its molecules and subject to the classical laws of 
thermodynamics(Davis, 1962; Shawagfeh, 1993). 



(16) 


subject to the boundary conditions 


y(0) = l, /(0) = 0, 


where m > 0 is constant. Substituting m - 0, 1,5 into Eq. (16) leads to the exact solution 




respectively. 

In other cases there aren’t any analytic exact solutions. Therefore, we apply ICSRBF method to 
solve the standard Lane-Emden Eq. (16) for m - 0, 1, 1.5, 2, 2.5, 3, 4 and 5. To this way now 
we can construct the residual functions as follows: 


Res(x) = xy®(^) + 2(y^^\x)) + x(y„(x)r. 


(17) 


As said before, to obtain the coefficient Res(x) is equalized to zero at N points by (11) : 


Res(xj) — 0, y = 0,1,2,..., A. 


By solving this set of equations, we can hnd the approximating function y„(x). 

Table 2 shows the comparison of the first zeros of standard Lane-Emden equations, from the 
present method and exact given by Horedt(Horedt, 2004) for m - 0, 1, 1.5, 2, 2.5, 3 and 4, 
respectively. 

Table 3, 4, 5 shows the approximation of y(x) for the standard Lane-Emden equation for m - 
2, 3, 4 respectively obtained by the method proposed in this paper and those obtained by 
Horedt(Horedt, 2004). 

Table 6 represents the ||Rei ||2 by the present method for m = 1.5, 2, 2.5, 3, 4 for several points 
and Table 7 represent the maximum error by the present method for m - 0, 1,5 for several 
points to how that the new method has appropriate convergence rate. The result graph of the 
standard Lane-Emden equation for m-0, 1, 1.5, 2, 2.5, 3, 4 and 5 is shown in Fig. 1. 

5.2. Example 2 (The isothermal gas sphere equation) 

For p(x) - 1, q(y(x)) - e^^^\ h(x) - 0, a - 2, A - 0 and B - 0, Eq.(13) is the isothermal 
gas sphere equation 



( 18 ) 




Figure 1: Graph of standai'd Lane-Emden equation for m=0, 1, 1.5, 2, 2.5, 3, 4 and 5. 


Table 8: Comparison of between present method and series solution given by Wazwaz(Wazwaz, 2001) and HFC 
method(Parand et al., 2010a) for isothermal gas sphere. 


X 

ICSRBF 

HFC 

Adomain 

0.00 

0.0000000000 

0.0000000000 

0.0000000000 

0.1 

-0.0016658338 

-0.0016664188 

-0.0016658339 

0.2 

-0.006533671 

-0.0066539713 

-0.0066533671 

0.5 

-0.0411539573 

-0.0411545150 

-0.0411539568 

1.0 

-0.1588276775 

-0.1588281737 

-0.1588273537 

1.5 

-0.3380194248 

-0.3380198308 

-0.3380131103 

2.0 

-0.5598230044 

-0.5598233120 

-0.5599626601 

2.5 

-0.8063408707 

-0.8063410846 

-0.8100196713 
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Figure 2: Graph of Isothermal gas sphere equation in comparison with Wazwaz solution(Wazwaz, 2001) 


subject to the boundary conditions 


y(0) = 0, /(0) = 0, 

Wazwaz(Wazwaz, 2001) by using ADM has produced a series solution as follow: 


, , 1 2 2 4 8 

yix) ^ —X H- X - 

6 5.4! 21.6! 


x® + 


122 

81.8!' 


61-67 ,0 
495.10! 


(19) 


We applied the ICSRBF method to solve this equation (18). We construct the residual function 
as follows: 


Res(x) - xyf\x) + 2(y^i^\x)) + (20) 

To obtain the coefficients ^, 5, Res{x) is equalized to zero at N point by (11) with y = 1.7: 

Res(xj)-0, j = 0,1,2, 

By solving this set of equations, we have the approximating function y(x). Table 8 shows the 
comparison of y(x) obtained by method proposed in this paper with (N = 40, r,,, = 6.5), and 
those obtained by Wazwaz(Wazwaz, 2001) (19) and HFC method by Parand et al(Parand et al., 
2010a). 

In order to compare the present method with those obtained by Wazwaz(Wazwaz, 2001) the 
resulting graph of Eq.(22) is shown in Fig. 2. The graph of residual for this equation for N - 
5, 10, 15, 20, 25, 30, 35 and 40 is shown in Fig. 3. 

5.3. Example 3 (The White Dwarf equation) 

For p{x) - 1, q(y(x)) - (y^ - cr)i, h{x) = 0, a - 2, A - I and B -0, Eq.(13) 
will be one of the white Dwarf equation that is absorbing to solve 

2 3 

y"{x) + -y'ix) + (y(x)^ - cr)^ = 0,x > 0, 

X 


( 21 ) 









Table 9: Comparison of y(:r), between present method and series solution given by Singh(Singh et al., 2009) and Haar 
method(Kaur et al., 2013) for White Dwarf equation. 


X 

(T 

ICSRBF 

MHAM 

Haar 

0.0001 

0.1 

0.999999 

0.999999 

1 

0.01 


0.999985 

0.999985 

0.999986 

0.1 


0.998579 

0.998578 

0.998581 

0.2 


0.994340 

0.994340 

0.994379 

0.4 


0.977738 

0.977738 

0.978345 

0.6 


0.951270 

0.951263 

0.953273 

0.7 


0.934823 

0.934801 

0.936497 

0.9 


0.896673 

0.896522 

0.895488 

0.0001 

0.2 

0.999999 

0.999999 

1 

0.01 


0.999988 

0.999988 

0.999988 

0.1 


0.998809 

0.998809 

0.998580 

0.2 


0.995255 

0.995251 

0.995296 

0.4 


0.981320 

0.981320 

0.981931 

0.6 


0.959049 

0.959045 

0.961128 

0.7 


0.945179 

0.945165 

0.947116 

0.9 


0.912913 

0.912812 

0.912833 

0.0001 

0.3 

0.999999 

0.999999 

1 

0.01 


0.999990 

0.999990 

0.999999 

0.1 


0.999025 

0.999025 

0.999028 

0.2 


0.996115 

0.996115 

0.996154 

0.4 


0.984690 

0.984690 

0.985295 

0.6 


0.966384 

0.966381 

0.968464 

0.7 


0.954953 

0.954944 

0.956942 

0.9 


0.928280 

0.928216 

0.928599 
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Figure 3: The ||/?C5||2 of Isothermal gas sphere equation in for N=5, 10, 15, 20, 25, 30, 35 and 40 


subject to the boundry conditions 


y(0) = l, y(0) = 0, 

which was introduced by Chandraskhar (Chandrasekhar, 1967) in his study of the gravitational 
potential of the degenerate White Dwarf stars. A series solution obtained by Singh(Singh et ah, 
2009) by using MHAM as follows: 

y{x) ^ 1 - - ;^trr^[5TiT^ + 14].*®. (22) 


where nr = Vl - cr. 

We applied the ICSRBF method to solve this equation (21). We construct the residual function 
as follows: 


Res(x) = xyf\x) + 2(y^^\x)) + x((y„(x)f - cr^)K (23) 

To obtain the coefficients ^iS, Res{x) is equalized to zero at N point by (11) with y - 1.5: 

Res(xj)-0, j = 0,1,2, 

By solving this set of equations,we have the approximating function y{x). Table 9 shows the 
comparison of y{x) obtained by method proposed in this paper with {N - 20, - 0.5), and 

those obtained by singh(Singh et al., 2009) (22) and Haar method by Kaur(Kaur et al., 2013). 
The result graph of the White Dwarf equation for cr- 0.1, 0.2 and 0.3 is shown in Fig. 4. 
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Figure 4: Graph of White Dwarf equation for cr=0. 1, 0.2 and 0.3 


Table 10: Comparison of between present method and series solution given by Wazwaz(Wazwaz, 2001) and HFC 
method(Parand et ah, 2010a) for Example 4. 


X 

ICSRBF 

HFC 

Adomain 

0.00 

1.0000000000 

1.0000000000 

1.0000000000 

0.1 

0.9980430038 

0.9981138095 

0.9980428414 

0.2 

0.9921896287 

0.9922758837 

0.9921894348 

0.5 

0.9519612468 

0.9520376245 

0.9519611019 

1.0 

0.8182430031 

0.8183047481 

0.8182516669 

1.5 

0.6254386159 

0.6254886192 

0.6258916077 

2.0 

0.4066218732 

0.4066479695 

0.4136691039 
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5.4. Example 4 

For p{x) - 1, q(y(x)) — sinh(y(x)), h{x) = 0, a - 2, A - I and B — 0, Eq.(13) will be one 
of the Lane-Emden type equations that is absorbing to solve 


y'{x) + -y{.x) + sinh(j) -0,x>0, 

X 


(24) 


subject to the boundry conditions 


y(0) = l, y(0) = 0, 

A series solution obtained by Wazwaz(Wazwaz, 2001) by using ADM is: 


^ ^ , (e^ - l)x^ 1 {e‘^ - \)x^ 

y(x) ^ 1-7^1-+ 


1 (2e® + 3e^ - Se'* - 2)x® 


30240 


1 (16e*^ - 104e® + 104^ - 61)x^ 


26127360 


(25) 


We applied the ICSRBE method to solve this equation (24). We construct the residual function 
as follows: 


Res(x) = xy^^x) + liy^^Kx)) + x sinh(y„(x)), (26) 

To obtain the coefficients Res{x) is equalized to zero at N point by (11) with y - 1.7: 

Res{xj)-0, j = 0,1,2, ...,N. 

By solving this set of equations,we have the approximating function y(x). Table 10 shows the 
comparison of y(x) obtained by method proposed in this paper with (N = 20, - l),and 

those obtained by Wazwaz(Wazwaz, 2001) (25) and HEC method by Parand et al.(Parand et al., 
2010a). 

In order to compare the present method with those obtained by Wazwaz(Wazwaz, 2001) the 
resulting graph of Eq.(22) is shown in Pig. 5. The graph of HResIb for this equation for N - 
5, 6, 7, 8, 9, 10, 13, 16 and 20 is shown in Fig. 6. 


5.5. Example 5 

For p(x) - 1, q(y(x)) - sin(y(x)), h(x) = 0, a = 2, A = 1 and B - 0, Eq.(13) will be one of 
the Lane-Emden type equations that is absorbing to solve 


y"(x) + -y'(x) -t- sin(y) = 0, x > 0, 

X 


( 27 ) 


subject to the boundry conditions 


y(0) = l, /(0) = 0, 
18 







Figure 5: Graph of equation of example 4 in comparing the present method and Wazwaz solution (Wazwaz, 2001). 
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Figure 6: The of example 4 for N = 5, 6, 7, 8, 9, 10, 13, 16 and 20 
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Table 11: Comparison of yCJc), between present method and series solution given by Wazwaz(Wazwaz, 2001) and HFC 
method(Parand et al., 2010a) for Example 5. 


X 

ICSRBE 

HEC 

Adomian 

0.0 

1.0000000000 

1.0000000000 

1.0000000000 

0.1 

0.9985979436 

0.9986051425 

0.9985979358 

0.2 

0.9943962892 

0.9944062706 

0.9943962733 

0.5 

0.9651778048 

0.9651881683 

0.9651777886 

1.0 

0.8636811302 

0.8636881301 

0.8636811027 

1.5 

0.7050451897 

0.7050524103 

0.7050419247 

2.0 

0.5064632371 

0.5064687568 

0.5063720330 


A series solution obtained by Wazwaz(Wazwaz, 2001) by using ADM is: 


y(x)- l-^kix^ + -l-kik 2 x'^ + ki( 

O 120 


3024 


rki.-- 


5040 


kl)x^ + kik2{- 


113 


3265920 




362880 


khx^ 


(28) 


where k\ - sin(l) and k 2 - cos(l). 

We applied the ICSRBF method to solve this equation (27). We construct the residual function 
as follows: 


Res{x) = xy^'nKx) + 2{y''^\x)) + x sin(y„(x)), (29) 

To obtain the coefficients Res{x) is equalized to zero at N point by (11) with y - 1.6: 

Res(xj)-Q, j = 0,1,2, 

By solving this set of equations, we have the approximating function y(ji:). Table 11 shows the 
comparison of y(x) obtained by method proposed in this paper with {N = 20, r^j - 2), and 
those obtained by Wazwaz(Wazwaz, 2001) (28) and HFC method by Parand et al. (Parand et al., 
2010a). 

In order to compare the present method with those obtained by Wazwaz(Wazwaz, 2001) the re¬ 
sulting graph of Eq.(27) is shown in Fig. 7. 

The graph of ||Rei ||2 for this equation for A = 5, 6, 7, 8, 9, 10, 13, 16 and 20 is shown in Fig. 8. 


6. Conclusion 

Lane-Emden equations occur in the theory of stellar structure and describe the temperature 
variation of a spherical gas cloud. White Dwarf equation appears in the gravitational potential of 
the degenerate White Dwarf stars. Lane-Emden type equations have been considered by many 
mathematicians as mentioned before(Dehghan and Shaken, 2008a). The fundamental goal of 
this paper has been to construct an approximation to the solution of non-linear Lane-Emden type 
equation in a semi-infinite interval. CSRBEs are proposed to provide an effective but simple way 
to improve the convergence of the solution by the collocation method. The validity of the method 
is based on the assumption that it converges by increasing the number of collocation points. A 
comparison is made among the numerical solution of Horedt(Horedt, 2004) and series solutions 
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Figure 7: Graph of equation of example 5 in comparing the present method and Wazwaz solution (Wazwaz, 2001). 
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Figure 8: The ||i?ej '||2 graph of example 5 for = 5, 6, 7, 8, 9, 10, 13, 16 and 20 
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of Wazwaz(Wazwaz, 2001), Singh et al.(Singh et al., 2009), Kaur et al.(Kaur et al., 2013) and the 
current work. It has been shown that our present work provides an acceptable approach for the 
Lane-Emden type equations. Also it was confirmed by maximum error and ||/?es ||2 figures, this 
approach has an exponentially convergence rate. Also, in this paper we show that the ICSRBF 
method for solving ordinary differential equations is simple and it has high accuracy and reliable 
convergence. 
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